Combining TP53 mutation and isoform has the potential to improve clinical practice.

The clinical importance of assessing and combining data on TP53 mutations and isoforms is discussed in this article. It gives a succinct overview of the structural makeup and key biological roles of the isoforms. It then provides a comprehensive summary of the roles that p53 isoforms play in cancer development, therapy response and resistance. The review provides a summary of studies demonstrating the role of p53 isoforms as potential prognostic indicators. It further provides evidence on how the presence of TP53 mutations may affect one or more of these activities and the association of p53 isoforms with clinicopathological data in various tumour types. The review gives insight into the present diagnostic hurdles for identifying TP53 isoforms and makes recommendations to improve their evaluation. In conclusion, this review offers suggestions for enhancing the identification and integration of TP53 isoforms in conjunction with mutation data within the clinical context.

A validated restriction enzyme ddPCR cg05575921 (AHRR) assay to accurately assess smoking exposure.

BACKGROUND & METHODS: In this study, a novel restriction enzyme (RE) digestion-based droplet digital polymerase chain reaction (ddPCR) assay was designed for cg005575921 within the AHRR gene body and compared with matching results obtained by bisulfite conversion (BIS) ddPCR and Illumina DNA methylation array. RESULTS: The RE ddPCR cg05575921 assay appeared concordant with BIS ddPCR (r2 = 0.94, P < 0.0001) and, when compared with the Illumina array, had significantly better smoking status classification performance for current versus never smoked (AUC 0.96 versus 0.93, P < 0.04) and current versus ex-smoker (AUC 0.88 versus 0.83, P < 0.04) comparisons. CONCLUSIONS: The RE ddPCR cg05575921 assay accurately predicts smoking status and could be a useful component of 'precision-medicine' chronic disease risk screening tools.

Biallelic ATP2B1 variants as a likely cause of a novel neurodevelopmental malformation syndrome with primary hypoparathyroidism.

ATP2B1 encodes plasma membrane calcium-transporting-ATPase1 and plays an essential role in maintaining intracellular calcium homeostasis that regulates diverse signaling pathways. Heterozygous de novo missense and truncating ATP2B1 variants are associated with a neurodevelopmental phenotype of variable expressivity. We describe a proband with distinctive craniofacial gestalt, Pierre-Robin sequence, neurodevelopmental and growth deficit, periventricular heterotopia, brachymesophalangy, cutaneous syndactyly, and persistent hypocalcemia from primary hypoparathyroidism. Proband-parent trio exome sequencing identified compound heterozygous ATP2B1 variants: a maternally inherited splice-site (c.3060+2 T > G) and paternally inherited missense c.2938 G > T; p.(Val980Leu). Reverse-transcription-PCR on the proband's fibroblast-derived mRNA showed aberrantly spliced ATP2B1 transcripts targeted for nonsense-mediated decay. All correctly-spliced ATP2B1 mRNA encoding p.(Val980Leu) functionally causes decreased cellular Ca2+ extrusion. Immunoblotting showed reduced fibroblast ATP2B1. We conclude that biallelic ATP2B1 variants are the likely cause of the proband's phenotype, strengthening the association of ATP2B1 as a neurodevelopmental gene and expanding the phenotypic characterization of a biallelic loss-of-function genotype.

Clonal dynamics limits detection of selection in tumour xenograft CRISPR/Cas9 screens.

Transplantable in vivo CRISPR/Cas9 knockout screens, in which cells are edited in vitro and inoculated into mice to form tumours, allow evaluation of gene function in a cancer model that incorporates the multicellular interactions of the tumour microenvironment. To improve our understanding of the key parameters for success with this method, we investigated the choice of cell line, mouse host, tumour harvesting timepoint and guide RNA (gRNA) library size. We found that high gRNA (80-95%) representation was maintained in a HCT116 subline transduced with the GeCKOv2 whole-genome gRNA library and transplanted into NSG mice when tumours were harvested at early (14 d) but not late time points (38-43 d). The decreased representation in older tumours was accompanied by large increases in variance in gRNA read counts, with notable expansion of a small number of random clones in each sample. The variable clonal dynamics resulted in a high level of 'noise' that limited the detection of gRNA-based selection. Using simulated datasets derived from our experimental data, we show that considerable reductions in count variance would be achieved with smaller library sizes. Based on our findings, we suggest a pathway to rationally design adequately powered in vivo CRISPR screens for successful evaluation of gene function.

Dynamic ctDNA Mutational Complexity in Patients with Melanoma Receiving Immunotherapy.

BACKGROUND: Circulating tumour DNA (ctDNA) analysis promises to improve the clinical care of people with cancer, address health inequities and guide translational research. This observational cohort study used ctDNA to follow 29 patients with advanced-stage cutaneous melanoma through multiple cycles of immunotherapy. METHOD: A melanoma-specific ctDNA next-generation sequencing (NGS) panel, droplet digital polymerase chain reaction (ddPCR) and mass spectrometry analysis were used to identify ctDNA mutations in longitudinal blood plasma samples from Aotearoa New Zealand (NZ) patients receiving immunotherapy for melanoma. These technologies were used in conjunction to identify the breadth and complexity of tumour genomic information that ctDNA analysis can reliably report. RESULTS: During the course of immunotherapy treatment, a high level of dynamic mutational complexity was identified in blood plasma, including multiple BRAF mutations in the same patient, clinically relevant BRAF mutations emerging through therapy and co-occurring sub-clonal BRAF and NRAS mutations. The technical validity of this ctDNA analysis was supported by high sample analysis-reanalysis concordance, as well as concordance between different ctDNA measurement technologies. In addition, we observed > 90% concordance in the detection of ctDNA when using cell-stabilising collection tubes followed by 7-day delayed processing, compared with standard EDTA blood collection protocols with rapid processing. We also found that the undetectability of ctDNA at a proportion of treatment cycles was associated with durable clinical benefit (DCB). CONCLUSION: We found that multiple ctDNA processing and analysis methods consistently identified complex longitudinal patterns of clinically relevant mutations, adding support for expanded clinical trials of this technology in a variety of oncology settings.

Complex Patterns of Genomic Heterogeneity Identified in 42 Tumor Samples and ctDNA of a Pulmonary Atypical Carcinoid Patient.

Tumor evolution underlies many challenges facing precision oncology, and improving our understanding has the potential to improve clinical care. This study represents a rare opportunity to study tumor heterogeneity and evolution in a patient with an understudied cancer type. A patient with pulmonary atypical carcinoid, a neuroendocrine tumor, metastatic to 90 sites, requested and consented to donate tissues for research. 42 tumor samples collected at rapid autopsy from 14 anatomically distinct sites were analyzed through DNA whole-exome sequencing and RNA sequencing, and five analyzed through linked-read sequencing. Targeted DNA sequencing was completed on two clinical tissue biopsies and one blood plasma sample. Chromosomal alterations and gene variants accumulated over time, and specific chromosomal alterations preceded the single predicted gene driver variant (ARID1A). At the time of autopsy, all sites shared the gain of one copy of Chr 5, loss of one copy of Chr 6 and 21, chromothripsis of one copy of Chr 11, and 39 small variants. Two tumor clones (carrying additional variants) were detected at metastatic sites, and occasionally in different regions of the same organ (e.g., within the pancreas). Circulating tumor DNA (ctDNA) sequencing detected shared tumor variants in the blood plasma and captured marked genomic heterogeneity, including all metastatic clones but few private tumor variants. This study describes genomic tumor evolution and dissemination of a pulmonary atypical carcinoid donated by a single generous patient. It highlights the critical role of chromosomal alterations in tumor initiation and explores the potential of ctDNA analysis to represent genomically heterogeneous disease. Significance: DNA sequencing data from tumor samples and blood plasma from a single patient highlighted the critical early role of chromosomal alterations in atypical carcinoid tumor development. Common tumor variants were readily detected in the blood plasma, unlike emerging tumor variants, which has implications for using ctDNA to capture cancer evolution.

Chromosomal Aberrations Accumulate during Metastasis of Virus-Negative Merkel Cell Carcinoma.

Merkel cell carcinoma is a rare, aggressive skin tumor initiated by polyomavirus integration or UV light DNA damage. In New Zealand, there is a propensity toward the UV-driven form (31 of 107, 29% virus positive). Using archival formalin-fixed, paraffin-embedded tissues, we report targeted DNA sequencing covering 246 cancer genes on 71 tumor tissues and 38 nonmalignant tissues from 37 individuals, with 33 of 37 being negative for the virus. Somatic variants of New Zealand virus-negative Merkel cell carcinomas partially overlapped with those reported overseas, including TP53 variants in all tumors and RB1, LRP1B, NOTCH1, and EPHA3/7 variants each found in over half of the cohort. Variants in genes not analyzed or reported in previous studies were also found. Cataloging variants in TP53 and RB1 from published datasets revealed a broad distribution across these genes. Chr 1p gain and Chr 3p loss were identified in around 50% of New Zealand virus-negative Merkel cell carcinomas, and RB1 loss of heterozygosity was found in 90% of cases. Copy number variants accumulate in most metastases. Virus-negative Merkel cell carcinomas have complex combinations of somatic DNA-sequence variants and copy number variants. They likely carry the small genomic changes permissive for metastasis from early tumor development; however, chromosomal alterations may contribute to driving metastatic progression.

Innate immune checkpoint inhibitor resistance is associated with melanoma sub-types exhibiting invasive and de-differentiated gene expression signatures.

Melanoma is a highly aggressive skin cancer, which, although highly immunogenic, frequently escapes the body's immune defences. Immune checkpoint inhibitors (ICI), such as anti-PD1, anti-PDL1, and anti-CTLA4 antibodies lead to reactivation of immune pathways, promoting rejection of melanoma. However, the benefits of ICI therapy remain limited to a relatively small proportion of patients who do not exhibit ICI resistance. Moreover, the precise mechanisms underlying innate and acquired ICI resistance remain unclear. Here, we have investigated differences in melanoma tissues in responder and non-responder patients to anti-PD1 therapy in terms of tumour and immune cell gene-associated signatures. We performed multi-omics investigations on melanoma tumour tissues, which were collected from patients before starting treatment with anti-PD1 immune checkpoint inhibitors. Patients were subsequently categorized into responders and non-responders to anti-PD1 therapy based on RECIST criteria. Multi-omics analyses included RNA-Seq and NanoString analysis. From RNA-Seq data we carried out HLA phenotyping as well as gene enrichment analysis, pathway enrichment analysis and immune cell deconvolution studies. Consistent with previous studies, our data showed that responders to anti-PD1 therapy had higher immune scores (median immune score for responders = 0.1335, median immune score for non-responders = 0.05426, p-value = 0.01, Mann-Whitney U two-tailed exact test) compared to the non-responders. Responder melanomas were more highly enriched with a combination of CD8+ T cells, dendritic cells (p-value = 0.03) and an M1 subtype of macrophages (p-value = 0.001). In addition, melanomas from responder patients exhibited a more differentiated gene expression pattern, with high proliferative- and low invasive-associated gene expression signatures, whereas tumours from non-responders exhibited high invasive- and frequently neural crest-like cell type gene expression signatures. Our findings suggest that non-responder melanomas to anti-PD1 therapy exhibit a de-differentiated gene expression signature, associated with poorer immune cell infiltration, which establishes a gene expression pattern characteristic of innate resistance to anti-PD1 therapy. Improved understanding of tumour-intrinsic gene expression patterns associated with response to anti-PD1 therapy will help to identify predictive biomarkers of ICI response and may help to identify new targets for anticancer treatment, especially with a capacity to function as adjuvants to improve ICI outcomes.

Impact of Clinical Data Veracity on Cancer Genomic Research.

Genomic analysis of tumors is transforming our understanding of cancer. However, although a great deal of attention is paid to the accuracy of the cancer genomic data itself, less attention has been paid to the accuracy of the associated clinical information that renders the genomic data useful for research. In this brief communication, we suggest that omissions and errors in clinical annotations have a major impact on the interpretation of cancer genomic data. We describe our discovery of annotation omissions and errors when reviewing an already carefully annotated colorectal cancer gene expression dataset from our laboratory. The potential importance of clinical annotation omissions and errors was then explored using simulation analyses with an independent genomic dataset. We suggest that the completeness and veracity of clinical annotations accompanying cancer genomic data require renewed focus by the oncology research community, when planning new collections and when interpreting existing cancer genomic data.

Survival of patients with small bowel neuroendocrine neoplasms in Auckland, Aotearoa New Zealand.

BACKGROUND: Small intestinal Neuroendocrine Neoplasms (SI-NENs) are the most common primary malignancy of the small bowel. The aim of this study is to define the survival of patients with an SI-NEN in Auckland, Aotearoa New Zealand (AoNZ). METHODS: A retrospective study of all patients diagnosed with a jejunal or ileal SI-NEN in the Auckland region between 2000 and 2012 was performed. The New Zealand NETwork! Registry was searched to identify the study cohort. Retrospective data collection was performed to collect stage, survival and follow up data. RESULTS: One hundred and seven patients were included in the study. The mean age of patients was 62.8 years (SD 11.9). The 5 and 10-year disease-specific survival for all patients was 66.1% (95% CI 56.5-75.7%) and 61.8% (95% CI 51.8-71.8%), respectively. Ten-year disease-specific survival was 100% for stage I and II, 74% (95%CI 61.7-84.4%) for stage III and 33.9% (95%CI 16.9-35.6%) for stage IV SI-NEN. Eleven of 40 (27.5%) patients with stage III disease had recurrence and 3 of 7 (42.8%) patients with stage IV disease had recurrence. In patients with stage IV disease, neither primary resection (HR 2.25, 95% CI 0.92-5.5) nor distant resection (HR 1.72, 95% CI 0.63-4.7) were significantly associated with a disease-specific or overall survival benefit. CONCLUSION: This study demonstrates that stage at SI-NEN diagnosis is associated with survival, but resection of the primary or distant metastases in patients with stage IV disease is not. There was no recurrence in patients with stage I or II disease after complete resection.

Tuberous sclerosis complex: a complex case.

Tuberous sclerosis complex (TSC) is an inheritable disorder characterized by the formation of benign yet disorganized tumors in multiple organ systems. Germline mutations in the TSC1 (hamartin) or more frequently TSC2 (tuberin) genes are causative for TSC. The malignant manifestations of TSC, pulmonary lymphangioleiomyomatosis (LAM) and renal angiomyolipoma (AML), may also occur as independent sporadic perivascular epithelial cell tumor (PEComa) characterized by somatic TSC2 mutations. Thus, discerning TSC from the copresentation of sporadic LAM and sporadic AML may be obscured in TSC patients lacking additional features. In this report, we present a case study on a single patient initially reported to have sporadic LAM and a mucinous duodenal adenocarcinoma deficient in DNA mismatch repair proteins. Moreover, the patient had a history of Wilms' tumor, which was reclassified as AML following the LAM diagnosis. Therefore, we investigated the origins and relatedness of these tumors. Using germline whole-genome sequencing, we identified a premature truncation in one of the patient's TSC2 alleles. Using immunohistochemistry, loss of tuberin expression was observed in AML and LAM tissue. However, no evidence of a somatic loss of heterozygosity or DNA methylation epimutations was observed at the TSC2 locus, suggesting alternate mechanisms may contribute to loss of the tumor suppressor protein. In the mucinous duodenal adenocarcinoma, no causative mutations were found in the DNA mismatch repair genes MLH1, MSH2, MSH6, or PMS2 Rather, clonal deconvolution analyses were used to identify mutations contributing to pathogenesis. This report highlights both the utility of using multiple sequencing techniques and the complexity of interpreting the data in a clinical context.

Clinically Relevant Biology of Hyaluronic Acid in the Desmoplastic Stroma of Pancreatic Ductal Adenocarcinoma.

ABSTRACT: Pancreatic ductal adenocarcinoma (PDAC) is notorious for its poor outcome. The presence of a dense desmoplastic stroma is a hallmark of this malignancy, and abundant hyaluronic acid (HA) within this stroma is a common feature of PDAC. At the end of 2019, an HA-targeting drug, after initial promise, failed phase 3 clinical trials in PDAC. This failure in the face of such strong evidence for biological importance forces us to turn back to the research and seek a better understanding of HA biology in PDAC. Therefore, in this review, we reexamine what is known about HA biology, the methods used to detect and quantify HA, and the ability of the biological models in which HA has been investigated to recapitulate an HA-rich desmoplastic tumor stroma. The role of HA in PDAC relies on its complex interplay with a range of HA-associated molecules, which have not been as extensively investigated as HA itself. Therefore, using large genomic data sets, we cataloged the abundance and activity in PDAC of molecules that modulate HA synthesis, degradation, protein interactions, and receptor binding. Based on their association with clinical characteristics and individual patient outcomes, we suggest a small number of HA-associated molecules that warrant further investigation as biomarkers and drug targets.

Bulk and Single-Cell Profiling of Breast Tumors Identifies TREM-1 as a Dominant Immune Suppressive Marker Associated With Poor Outcomes.

BACKGROUND: Triggering receptor expressed on myeloid cells (TREM)-1 is a key mediator of innate immunity previously associated with the severity of inflammatory disorders, and more recently, the inferior survival of lung and liver cancer patients. Here, we investigated the prognostic impact and immunological correlates of TREM1 expression in breast tumors. METHODS: Breast tumor microarray and RNAseq expression profiles (n=4,364 tumors) were analyzed for associations between gene expression, tumor immune subtypes, distant metastasis-free survival (DMFS) and clinical response to neoadjuvant chemotherapy (NAC). Single-cell (sc)RNAseq was performed using the 10X Genomics platform. Statistical associations were assessed by logistic regression, Cox regression, Kaplan-Meier analysis, Spearman correlation, Student's t-test and Chi-square test. RESULTS: In pre-treatment biopsies, TREM1 and known TREM-1 inducible cytokines (IL1B, IL8) were discovered by a statistical ranking procedure as top genes for which high expression was associated with reduced response to NAC, but only in the context of immunologically "hot" tumors otherwise associated with a high NAC response rate. In surgical specimens, TREM1 expression varied among tumor molecular subtypes, with highest expression in the more aggressive subtypes (Basal-like, HER2-E). High TREM1 significantly and reproducibly associated with inferior distant metastasis-free survival (DMFS), independent of conventional prognostic markers. Notably, the association between high TREM1 and inferior DMFS was most prominent in the subset of immunogenic tumors that exhibited the immunologically hot phenotype and otherwise associated with superior DMFS. Further observations from bulk and single-cell RNAseq analyses indicated that TREM1 expression was significantly enriched in polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs) and M2-like macrophages, and correlated with downstream transcriptional targets of TREM-1 (IL8, IL-1B, IL6, MCP-1, SPP1, IL1RN, INHBA) which have been previously associated with pro-tumorigenic and immunosuppressive functions. CONCLUSIONS: Together, these findings indicate that increased TREM1 expression is prognostic of inferior breast cancer outcomes and may contribute to myeloid-mediated breast cancer progression and immune suppression.

Breast Cancer Patient Prognosis Is Determined by the Interplay between TP53 Mutation and Alternative Transcript Expression: Insights from TP53 Long Amplicon Digital PCR Assays.

The TP53 gene locus is capable of producing multiple RNA transcripts encoding the different p53 protein isoforms. We recently described multiplex long amplicon droplet digital PCR (ddPCR) assays to quantify seven of eight TP53 reference transcripts in human tumors. Here, we describe a new long amplicon ddPCR assay to quantify expression of the eighth TP53 reference transcript encoding ∆40p53α. We then applied these assays, alongside DNA sequencing of the TP53 gene locus, to tumors from a cohort of New Zealand (NZ) breast cancer patients. We found a high prevalence of mutations at TP53 splice sites in the NZ breast cancer cohort. Mutations at TP53 intron 4 splice sites were associated with overexpression of ∆133TP53 transcripts. Cox proportional hazards survival analysis showed that interplay between TP53 mutation status and expression of TP53 transcript variants was significantly associated with patient outcome, over and above standard clinical and pathological information. In particular, patients with no TP53 mutation and a low ratio of TP53 transcripts t2 to t1, which derive from alternative intron 1 acceptor splice sites, had a remarkably good outcome. We suggest that this type of analysis, integrating mutation and transcript expression, provides a step-change in our understanding of TP53 in cancer.

Genomic and signalling pathway characterization of the NZM panel of melanoma cell lines: A valuable model for studying the impact of genetic diversity in melanoma.

Melanoma is a disease associated with a very high mutation burden and thus the possibility of a diverse range of oncogenic mechanisms that allow it to evade therapeutic interventions and the immune system. Here, we describe the characterization of a panel of 102 cell lines from metastatic melanomas (the NZM lines), including using whole-exome and RNA sequencing to analyse genetic variants and gene expression changes in a subset of this panel. Lines possessing all major melanoma genotypes were identified, and hierarchical clustering of gene expression profiles revealed four broad subgroups of cell lines. Immunogenotyping identified a range of HLA haplotypes as well as expression of neoantigens and cancer-testis antigens in the lines. Together, these characteristics make the NZM panel a valuable resource for cell-based, immunological and xenograft studies to better understand the diversity of melanoma biology and the responses of melanoma to therapeutic interventions.

Circulating tumor DNA is a sensitive marker for routine monitoring of treatment response in advanced colorectal cancer.

Accurate assessment of chemotherapy response provides the means to terminate ineffective treatment, trial alternative drug regimens or schedules and reduce dose to minimize toxicity. Here, we have compared circulating tumor DNA (ctDNA) with carcinoembryonic antigen (CEA) for the cycle by cycle assessment of chemotherapy response in 30 patients with metastatic colorectal cancer. CtDNA (quantified using individualized digital droplet PCR (ddPCR) assays) and CEA levels were determined immediately prior to each chemotherapy cycle over time periods ranging from 42-548 days (average of 10 time points/patient). Twenty-nine/thirty (97%) patients had detectable ctDNA compared with 83% whose tumors were CEA-positive (>5 ng/ml) during the monitoring course. Over the course of treatment, 20 disease progression events were detected by computed tomography; ctDNA predicted significantly more of these events than CEA (16 (80%) versus 6 (30%), respectively; P-value = 0.004). When progression was detected by both ctDNA and CEA, the rise in ctDNA occurred significantly earlier than CEA (P-value = 0.046). Partial responses to chemotherapy were also detected more frequently by ctDNA, although this was not significant (P-value = 0.07). In addition, another 28 colorectal cancer patients who underwent potentially curative surgery and showed no evidence of residual disease were monitored with ctDNA for up to 2 years. Clinical relapse was observed in 6/28 (21%) patients. Four out of 6 of these patients showed a significant increase in ctDNA at or prior to relapse. Overall, ctDNA analyses were able to be performed in a clinically relevant timeline and were a more sensitive and responsive measure of tumor burden than CEA.

Specialized Cell-Free DNA Blood Collection Tubes Can Be Repurposed for Extracellular Vesicle Isolation: A Pilot Study.

Background: Liquid biopsies offer a minimally invasive approach to patient disease diagnosis and monitoring. However, these are highly affected by preprocessing variables with many protocols designed for downstream analysis of a single molecular biomarker. Here we investigate whether specialized blood tubes could be repurposed for the analysis of an increasingly valuable biomarker, extracellular vesicles (EVs). Methods: Blood was collected from three donors into K3-EDTA, Roche, or Streck cell-free DNA (cfDNA) collection tubes and processed using sequential centrifugation either immediately or after storage for 3 days. MicroEV were collected from platelet-poor plasma by 10,000 g centrifugation and NanoEVs isolated using size exclusion chromatography. Particle size and counts were assessed by Nanoparticle Tracking Analysis, protein quantitation by bicinchoninic acid assay (BCA) assay, and dot blotting for blood cell surface proteins. Results: MicroEVs and NanoEVs could be isolated from plasma collected using all three tube types. Major variations were seen with delayed time to processing. Both MicroEV particle number and protein content increased with the processing delay. The NanoEV number did not change with the time-delay but their protein quantity increased. EV-associated proteins predominantly arose from platelets (CD61) and erythrocytes (CD235a). However, leukocyte marker CD45 was only increased in NanoEVs from ethylenediaminetetraacetic acid (EDTA) tubes, suggestive of stabilization of nucleated cells by the specialized blood tubes. Epithelial cell surface marker EpCAM, often used as a marker of cancer, remained the same across conditions in both MicroEV and NanoEV preparations indicating that these EVs were stable with time. Conclusions: Specialized cfDNA collection tubes can be repurposed for MicroEV and NanoEV analysis; however, simple counting or using protein quantity as a surrogate of EV number may be confounded by preanalytical processing. The EVs would be suitable for disease selective EV subtype analysis if the molecular target of interest is not present in blood cells.

PDX1 DNA Methylation Distinguishes Two Subtypes of Pancreatic Neuroendocrine Neoplasms with a Different Prognosis.

DNA methylation is a crucial epigenetic mechanism for gene expression regulation and cell differentiation. Furthermore, it was found to play a major role in multiple pathological processes, including cancer. In pancreatic neuroendocrine neoplasms (PNENs), epigenetic deregulation is also considered to be of significance, as the most frequently mutated genes have an important function in epigenetic regulation. However, the exact changes in DNA methylation between PNENs and the endocrine cells of the pancreas, their likely cell-of-origin, remain largely unknown. Recently, two subtypes of PNENs have been described which were linked to cell-of-origin and have a different prognosis. A difference in the expression of the transcription factor PDX1 was one of the key molecular differences. In this study, we performed an exploratory genome-wide DNA methylation analysis using Infinium Methylation EPIC arrays (Illumina) on 26 PNENs and pancreatic islets of five healthy donors. In addition, the methylation profile of the PDX1 region was used to perform subtyping in a global cohort of 83 PNEN, 2 healthy alpha cell and 3 healthy beta cell samples. In our exploratory analysis, we identified 26,759 differentially methylated CpGs and 79 differentially methylated regions. The gene set enrichment analysis highlighted several interesting pathways targeted by altered DNA methylation, including MAPK, platelet-related and immune system-related pathways. Using the PDX1 methylation in 83 PNEN, 2 healthy alpha cell and 3 healthy beta cell samples, two subtypes were identified, subtypes A and B, which were similar to alpha and beta cells, respectively. These subtypes had different clinicopathological characteristics, a different pattern of chromosomal alterations and a different prognosis, with subtype A having a significantly worse prognosis compared with subtype B (HR 0.22 [95% CI: 0.051-0.95], p = 0.043). Hence, this study demonstrates that several cancer-related pathways are differently methylated between PNENs and normal islet cells. In addition, we validated the use of the PDX1 methylation status for the subtyping of PNENs and its prognostic importance.

An estimate of limited duration cancer prevalence in New Zealand using 'big' data.

AIMS: Increases in cancer survival may increase cancer prevalence and demand for healthcare. We aimed to estimate cancer prevalence in the New Zealand population. METHODS: We used national linked health, social and census datasets from the Stats NZ Integrated Data Infrastructure to identify the number of New Zealand residents who had at least one cancer diagnosis in New Zealand. We included all primary cancers recorded on the New Zealand Cancer Registry from January 1995 to June 2013, and used the 2013 census for demographic and socioeconomic data. RESULTS: On 30 June 2013, 140,600 of 4,438,900 (3.2%) New Zealand residents had been diagnosed with cancer in the last 18.5 years. In ≥15 year olds, the age-standardised prevalence of cancer diagnosed 0 to ≤1 year, and >1 to ≤5 years, prior to 30 June 2013 was 0.4% and 1.1% in men and 0.3% and 0.9% in women, respectively. Over the 18.5-year period prevalence was greatest in the oldest ages, European/Other, highest qualified, highest income, least deprived, ex-smokers, and Canterbury, Bay of Plenty and Nelson/Marlborough District Health Boards (age-standardised). CONCLUSIONS: Groups with the highest survival and the greatest access to healthcare had the highest cancer prevalences.